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Abstract 

In this paper we propose a numerical scheme for the class of backward doubly 
stochastic (BDSDEs) with possible path-dependent terminal values. We prove that 
our scheme converge in the strong L 2 -sense and derive its rate of convergence. As 
an intermediate step we derive an L 2 -type regularity of the solution to such BDSDEs. 
Such a notion of regularity which can be though of as the modulus of continuity of the 
paths in an L 2 -sense, is new. 
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1 Introduction 

In this paper we are interested in the following backward doubly stochastic differential 
equations (BDSDEs, in short): 

Y t = £ + f f{s,X s ,Y s ,Z s )ds+ f T g{s,X s ,Y s )d~B s - J* Z s dW s , (1.1) 

where W and B is two independent Brownian motion defined on (Q, f ,P) a complete prob- 
ability space. This kind of equation has two different directions of stochastic integrals: stan- 
dard (forward) stochastic integral driven by W and backward stochastic one driven by B. 
Initiated by Pardoux-Peng ||20l , BDSDEs is connected to quasi-linear stochastic partial dif- 
ferential equations (SPDEs, in short) in order to derive Feynman-Kac formula for SPDEs. 
In this setting, BDSDEs have been extensively studied in the past decade. We refer the read- 
ers to the papers of Buckdahn and Ma 00 [H, Aman et al. [2], Aman [ 1 1, Bahlali et al.f3l. 
and reference therein for more information on both theory and applications, especially in 
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mathematical finance and stochastic control, for such equations. In contrast, there was little 
progress made in the direction of the numerical implementation of BDSDEs. In special 
case of BDSDEs (g = 0) called BSDEs, many efforts have been made in this direction as 
well. 

Up to now basically two types of schemes have been considered. Based on the theo- 
retical four-step scheme from |[T8l . first type of numerical algorithms for BSDEs have been 
developed by Douglas et al. [11] and more recently by Milstein and Tretyakov [19]. The 
main focus of these algorithms is the numerical solution of parabolic PDEs which is related 
to BSDEs. 

A second type of algorithms works backwards through time and tries to tackle the 
stochastic problem directly. Bally [4] and Chevance [10] were the first to study this type 
of algorithm with a (hardly implementable) random time partition under strong regularity 
assumptions. The works of Ma et al. [ 16 ] and Briand et al. Q ares in the same spirit, re- 
placing, however, the Brownian motion by a binary random walk. Recently, Zhang proved, 
in ETI . a new notion of //-regularity on the control part Z of the solution which allowed 
proof of convergence of this backward approach with deterministic partitions under rather 
weak regularity assumptions (see [El J, [6 ] and [ 13 ]) for different algorithms. All numerical 
schemes provide alternative to construct algorithm for PDEs. To the best of our knowledge 
to date there has been no discussion in the literature concerning numerical algorithms in 
the spirit of the last three works cited above in the general case i.e g ^ 0. This constitutes 
an insufficiency when we know that almost all the deterministic problems in these applied 
fields (PDEs) have their stochastic counterparts (SPDEs). 

In this paper, to correct this empty, our goal is to build a numerical scheme following 
the idea used by Bouchard and Touzi [6] and study its convergence. These results were 
important from a pure mathematical point of view as well as from application in the world 
of finance. Particulary, this numerical scheme opens the way for possible algorithm for 
determining the price of options on financial assets whose dynamics is solution of SPDEs. 

Similarly to the special case g = 0, the main difficulty lies in the approximation of 
the "martingale integrand" Z. In fact, in a sense the problem often comes down to the 
path regularity of Z. However, in case g ^ 0, this regularity becomes a natural question 
to ask. Therefore, the first main result in this paper is to derive the path regularity called 
L 2 -regularity for BDSDEs with the terminal value £, is the form <&(X), where X and <!>(.) 
are respectively diffusion process and L°°-Lipschitz functional (see Section 3 for precise 
definition). The proof is heavily related to Girsanov's Transformation which exists in the 
BDSDEs case only if g do not depend to z. 

The above L 2 -regularity result allow us to provide the rate of convergence of our numer- 
ical scheme which is different from the one constructed in O. Indeed, since BDSDEs have 
two directions of integral, our numerical scheme need at each step the conditional expecta- 
tion with respect the filtration defined by = c(X™, j < i)Vc(B tj , j < i). However 
we obtain the same convergence rate. 

The rest of this paper is organized as follows. In Section 2, we introduce some fun- 
damental knowledge and assumptions of BDSDEs. Section 3 is devoted to L 2 -regularity 
results. In Section 4, we built our numerical scheme and prove the rate of convergence. 
Finally in section 5, we focus some ideas for the regression approximation and give it con- 
vergence rate. 
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2 Preliminaries 



Let (Q, f ,P) be a complete probability spaces, and T > be fixed throughout this paper. 
Let {W t ,0 <t<T} and {fi ? ,0 < t < T} be two mutually independent standard Brownian 
motions processes, with values respectively in R d and R e , defined on (£2,!F ,P). Let 5\£ 
denote the class of P-null sets of J . For each t G [0, T], we define 



■ w . 



where for any process (r| v : < 5 < T) , fF^ = o{r\ r — r\ s ,s<r<t}\/c\C, 7^ = 7§ v 

We note that since the collection {f t )t>o is neither increasing nor decreasing, it does 
not constitute a filtration. Therefore, we define the filtration (fF t ) t >o by 



ft 



which contains f t and play a key role in the build of our numerical scheme. 

For any real p > 2 and k G N*, let S p (M. k ) denote the set of jointly measurable processes 
{Xt}te{o,T} taking values in R k which satisfy 



(0 11*1 



E 



sup \X t \ p 

0<t<T 



< +°°; 



(ii) X t is jF r measurable, for any t G [0, T]. 

We denote similarly by M P (R ) the set of (classes of d¥®dt a.e. equal) ^-dimensional 
jointly measurable processes which satisfy 

i 

\ 9 

|2, 



(011*1 



MP 



dt 



< +°°; 



R, such that for some constant 



(»)*? is ir f -measurable, for a.e. t G [0,T]. 
We denote by 

• W 1,00 (R*) the space of all measurable functions \\f : R k — > 

^>0 it holds that \y(x) -\\r(y)\ <K\x-y\, Vx,vG 

• B the space of all cadlag functions defined on [0, T\ ; 

• C™([0, T] x R k ) the space of all continuous functions (not necessary bounded) \|/ : [0, T] x 

R k — > R, such that \|/ has uniformly bounded derivatives with respect to the spatial 
variables up to order m. We often denote C™ = C™([0, T\ x R fc ) for simplicity, when 
the context is clear. 



Let 



b 
a 

f 
8 



[0,71 x 
[0,T]x 
[0,T}x 
[0,T] xR d x\ 



l d , 
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X K X . 
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be the functions satisfying the following assumptions: there exists constant K > such that 

for all s,s' G [0,T],x,x' G R d , y,y' G R, z,z' G R d , 

(HI) |fc(j,jc) - b(s,x')\ + ||a(j,x) - a(s,x')\\ < K\x-x'\. 
(H2) 

(0 \f{s,x } y,z)- f(s\x\y\z^\ 2 <K{\s-sf +\x-xf + \y-yf +\z-zf) , 

(ii) \g(s,x,y)-g(s,x',y')\ 2 <K(\s-s'\ 2 + \x-x'\ 2 +\y-y'\ 2 ). 
(H3) sup < f < r {|^,0)| + |o(f,0)| + 1^,0,0,0)1 + |*(f,0,0)|} < K. 
Given £ G L 2 (£2, jF r ,P;M c/ )» denote (X,F,Z) be the solution to the following FBDSDE: 

X, = x+ f b(s,X s )ds + f a(s,X s )dW s (2.1) 
Jo Jo 

Y t = + J f(s,X s ,Y s ,Z s )ds + J g(s,X s ,Y s )dT s - j Z s dW s . (2.2) 

Let now recall some standard results appear in SDEs and BDSDEs literature. 

Proposition 2.1 (Karatzas and Shreve |[T4l ). Assume (HI) holds. Then for any initial con- 
dition x G M. d , FSDE (12.11 ) has a unique solution (X t )o< t <T belong to S p (M. d ). 

Moreover, for any p > 2, there exists a constant C p > 0, depending only on T, K and p, 
such that 

E y sup \X t \Pj <C p (^\x\p + £[\b(t,0)\P + \a(t,0)\ p ]dt^j 



and 



E[\X t -X s \ p ] <C p E[ \x\ p + sup \b(t,0)\ p + sup \a(t,0)\ p ) \t-s 

0<t<T 0<t<T 



P/2 



Proposition 2.2 (Pardoux and Peng EDI ). Under assumption (H2), BDSDE \2.2\ has a 

unique solution (Y t ,Z t )o< t <T in 5 P (M) x M p {W i ). 

Moreover, for any p > 1, there exists a constant C p > 0, depending only on T, K and p, 
such that 



E sup \Y t 

\o<f<r 



and 



E[\Y t -Y s \P] <C P E- 



' t \ p + ^ \Z s \ 2 dsy 2 ^j <C p E^\ p + £ '[|/(f,0,0,0)r + \g(t,0,0)\ p ]dt 



\t,\ p + sup \f(t,0,0,0)\ p + sup \g(t,0,0)\ p 

0<t<T 0<t<T 



\t-s\ p - 1 + 



\ p/2" 



Remark 2.3. In Proposition 12. 2[ existence and uniqueness result needs only Lipschitz con- 
dition on / and g with respect variables y and z uniformly in t and x. 
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Proposition 2.4 (Stability). Let (Z E ,7 e ,Z E ) be the solution to the perturbed FBDSDE (l2~TT ) 
and (12.21 ) in which the coefficients are replaced by b E ,G £ , f E , g E , with initial state x E an 
terminal value Assume that the assumption (HI) and (H2) hold for all coefficients 
b e , <T E , / E , g E and assume that lim E ^o* e = x , and for fixed (x,y,z), 

limEj [ [\b e (t,x)-b(t,x)\ 2 + \a e (t,x) -a(t,x)\ 2 ]dt\ = 0, 

e^O [Jq J 
limE||^-^| 2 + j o T [\f(t,x,y,z)-f(t,x,y,z)\ 2 + \g £ (t,x,y)-g(t,x,y)\ 2 ]dt} = 0. 

Then, we have 

lim E i sup |X* -X,\ 2 + sup |F, E - Y t \ 2 + sup |Zf - Z t \ 2 \ = 0. 
lo<f<r o<r<r o<f<r J 

3 L 2 -regularity result for BDSDEs 

In this section we establish the first main result of this paper which is L 2 -regularity of the 
martingale integrand Z; and can be thought of as the modulus of continuity of the paths in 
an L 2 sense. Such a regularity, combined with the estimate for X and Y, plays a key role for 
deriving the rate of convergence of our numerical scheme in Section 4. We shall consider a 
class of BDSDEs with terminal values which are path dependent i.e of the form ^ = <&(X), 
where a deterministic functional <J> : D — > R satisfies: 

(H4) (L°°-Lipschitz condition), there exists a constant K such that 

\^(X r )-^(X 2 )\<K sup |Xi(0-X 2 (t)|, VXi,X 2 GD. (3.1) 

o<f<r 

(H5) ^(0) is bounded by K, where denotes the constant function taking value on [0, T]. 

This approximation due to Ma and Zhang lfT7l . for L°°-Lipschitz functional will be useful 
in the sequel. 

Lemma 3.1. Suppose (H4) and (H5) hold. Let Yl = {n} be a family of partitions of [0, T]. 
Then there exists a family of discrete functionals {h % : 71 € IT} such that 

(i) for each 7t G n, assuming K : = to < t\ < t 2 < ■ ■ ■ < t n = T, we have that h n £ 
c i( R d(*+i)^ anc i satisfies: 

j^\d Xi h % (x)\<K, y,x=(x ,x u --;x n )eR d{n+1 \ (3.2) 
i=\ 

where K is the same constant as that in (13.11) . 
(h) for any XgD, it holds that 

lim \h n (X to ,X tl ,--;X tn )-<S>(X)\=0, (3.3) 

where \%\ = maxi<,-<„ If; — 
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Our main result in this section is the following theorem. 

Theorem 3.2. Assume (H1)-(H5). Let Kq : sq < ■ ■ -,s m be any partition of [0,T], and for 
each 1 < i < m, let define 



i 



Si-Si-l 



-E 



Z s ds\j ; 



Si- I 



S|-l 



Then there exists a constant C depending only on T and K, such that 



E 



max sup \Y t -Y s . f + Y / \Z S -Zf \ 2 ds 

l<i<'% M <K Sl ' i?\JSi-i " 1 



<C|tuo| 



(3.4) 



(3.5) 



In the sequel, let % : = to, ■ • ■, t n = T be any partition of [0, T] finer than 7to and without 
loss of generality, we assume Si = t/ j for i = l,...,m. Since <!> satisfies the L°°-Lipschitz 
condition (13.11 ). by virtue of Lemma I3TT1 one can find h % G C 1 (R rf(n+1 ') satisfying (13.21 ) and 
(T33T) . Let (F^Z 71 ) be the solution to BDSDE: 

F/ 1 = h^(X t0 ,...,X t J+J t f( S ,X s ^,Z^)d S + J t g( S ,X s J^d¥ s - j Z*dW,{3.6) 

Moreover, setting S n = (Z n ,Z % ), with E n = {X,Y n ), let (VX, VT 71 , V'Z 71 ) be the unique 
solution of the following FBDSDE: 

VXt = I d + f b x {r,X r )VXrdr+ f G x (r,X r )VX r dW r , 
Jo Jo 

n 7)h % r T 
V'F* = 1 £—(X t0 ,...,X tn )VX tj + lf x (®?)VX r + fii®«)V'Y? + U&*)V'Z?}dr (3.7) 

+ j\g x {Z n r )VX r + gy {Z n r )VX)\d~B r - J^VZfdWr, te[t u T], i = 0,...,n-l. 



We denote 



where 



^° = [ T M®«)VX r N; 1 dr+ f 1 g x (Z*)VX r N; l tB r ; 
Jo Jo 

= h K (X t0 ,---,X tn )VX ti N T \i = \,--;n, 

N t = exp^f y (e^dr + ^g y (E^ r -y^\g y (E^\ 2 dry 

M t = e xpU*f z (<^)dW r -y^\M®?)\ 2 dr\. (3.8) 

The following technical lemma is the building block of the proof of Theorem |3.2j 

Lemma 3.3. Let consider the partition ji defined above and h n given by Lemma IJ.il and 
assume G, b, f, g, E C\. Then for all i = 1, ...,n 

V% n = U° + £§A Mr 1 ^ - £f x (e*)VX r N- l drN t - f\ x {Z*)VX r N- l d~B r N u 
where % = E (M r ^>/) , j = 0,---,n. 
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Proof. For each < i < n, we recall (VX, VT 71 , V'Z 71 ), the solution of the linear FBDSDE 
(I3J1) . Let (fX°) and (Y 7 ,^), J = 1, ■ ' be the solution of the BDSDEs 



it = J^fM)VXr + fy(&M + fM)C]dr 

+ f\g x {Z*)VX r + g y (^)f]tB r - £ $dW r ; (3.9) 
it = ^(X t0 , ,X tl ,)VX tj +J t [f y {®*)t + fzm%l\dr 

respectively, then we have the following decomposition: 

V ! 'F/ = 7? + L^, Je[M4 (3-10) 

/=« 

Recall (13.81 ) and since and g y are uniformly bounded, by Girsanov's Theorem (see, 
e.g., HI) weknowthatMisaP-martingaleon [0,T], and W t = W t — j' f z (®f)dr, t G [0,T], 
is an J,-Brownian motion on the new probability space {£l,!F ,P), where P is defined by 

d¥ mT - 

Now for <i <n, define 

% = i t N-\ ll = QN t -\ te[0,T}. 
Then, using integration by parts and equation (13.91 ) we have 

Therefore, by the Bayes rule (see e.g, lfl4l Lemma 3.5.3) we have, for t G [0, T] 

i=%N t = E(M T t, i \f t )M- l N t =^M- l N t , 
where for i = 0, • • -,n, 

% = E{M T tj\? t } = E(M T ?) + [\\dW s . (3.11) 

Jo 

Recalling (I3T23T) . M T G LP(f T ) and VX G Z/(F) for all p > 2. 
Therefore for each p > 1 . (13.21 ) leads to 

e(v |AfrVl} <C p e{\M t \ p sup |VX t | p l. (3.12) 

U=i J i °< ? < r J 

In particular, for each j, MjlJ G L 2 (!Ft), thus (13.1 II ) makes sense. Finally the result follows 
by (T3T01) . □ 



7 



Proof of Theorem [O] Let us recall the partition 7to defined above and consider \tzq\ it mesh 
defined by 



7to = max \Si — Si-i\. 

0<i<m 



Applying Proposition 12.21 we get 

E(\Y t -Y Si _f) <C\n \, t€[si-i,Si), i = l,---,m, 
which together with Burkolder-Davis-Gundy inequality implies 



E 



max sup \Y t — Y Sil 
i<i<m Sjl < t < Si 



<C\k \. 



(3.13) 



The estimate for Z is a little involved. This part will be divide in two steps. 
Step 1 

First we assume that b, G, f,g£ C\. It follows from Lemma |3~T1 together with Proposition 
2.4 that 



lim e( sup \Y t n -Y t \ 2 + f \Zf -Z t \ z dt\ = 0. 
|to|->o U<f<r Jo J 



(3.14) 



On the other hand, according to CO}, Z£° , € L 2 (D., f s ._ { ). Then since Z£_, G L 2 (H, ), 
it follows from Lemma 3.4.2 of 11221 . page 71, that 



E 



< 2E 



£/ |z,-2* x | 2 <fa 



< E 



1=1^-1 



Z<-Z*J 2 ds 



£/ (|z J -z-| 2 +|z?-z^_ 1 | 2 )^ 



< C|7I |+E 
By (I3.14I) and (I3.15I ). it remains to prove that 



Z7-Z5 fds 



(3.15) 



Z?-Z? \ 2 ds 

.) 1 I 



<C|7I |, 



(3.16) 



where C is independent of 71 or 7to. Now we fix z'q. For f G [jj _i,j ; - ), it follows from 
Proposition 2.3 of [20] together with Lemma 3.3 that 



Zf 



# + I#J M- 1 - Jj x (e^)VX r N- l dr- £ g^VX.N-'dt,- 



N t [VX t ]- l a(X t 



Therefore, 



|Z*-Z* i |<//+/ 2 + /3+^ 



(3.17) 
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where (recalling that j,- _i = t\ j 



M7 1 JV,. ,[VX S . .TWx,. 



M-%[VX t ]- l c(X t ) -M-^N^ VX,, : ,] l a(X s . ^ 



/, 4 = Kl. 



with 



Ai 



MQ^VXrN-'dr N t [VX t ]- l a(X t ) 



and 



jf &(E?) VX^K) Ai [VX] - 1 a(X f ) 

g^S^VX.iV- 1 ^ AT S . [VX s , o _ 1 ]- 1 a(X., o _ 1 



Recalling (13 - 8b - and noting that f y , f z and g y are uniformly bounded, one can deduce 
that, for all p > 1, there exists a constant C p depending only on T, K and p, such that 



sup W + IAT'M < C p ; E( sup [\M t \ p + \M; l \ p ] ) <C P ; 
Kf<r / Vo<f<r / 

E ( I - A^v | p + | A^r 1 - 1 1 p ) < C^-s^ 2 ; 



(3.18) 



E(|M f -M J | p + |M f _1 -M7 1 | p ) < C p \t-s\ p l 2 . 
Thus, applying Proposition 2.1 and 2.2 one can show that 

E(|/, 3 | 2 )<C|7i |, 

and 

E(|/, 4 | 2 )<C|7io|, 
Recalling (13.111 ) and (13.21 ) we have 

\% + T%\<Ce\ sup WX t \f t 
j>i lo<f<r 



(3.19) 



(3.20) 



Thus by using again Proposition 2.1 and 2.2 together with (13.231) . we get 

E(\lf\ 2 ) <C[jto|. 
As proved in ET1 (see proof of theorem 3.1), we have 

E(|//| 2 )<C|7i |. 



(3.21) 



(3.22) 



Combining (l3~T9l) . (f3T20b . (f3T2Tb and (13^221 . we deduce from (|3~T7T ) that (l3~T6l) holds, which 
ends the proof for the smooth case. 
Step 2 

Let consider the general case i.e b, a, /, g are only Lipschitz. For cp = b, a, f, g, it not 
difficult to construct via a convolution method, for any £ > 0, the function cp E € C\ be the 
smooth mollifiers of cp such that the derivatives of cp E are uniformly bounded by K and 
limg^oq) 8 = cp. Let (X E ,Y E ,Z E ) and (X E ,Y n > E ,Z 7C,E ) denote the solution to corresponding 
FBDSDE replaced cp by cp e and set 

N E = exp(^ e (0: t - £ Mr + ^ 

Mf = e X p|£^(0^)^ r -i^^(0^)| 2 Jr|. 

Then one can derive, since the function f£ ,/ z E and g E are uniformly bounded by K, with the 
standard calculus about BSDEs, that, for all p > 1, there exists a constant C p independent 
on £ (depending only on T, K and p), such that 

e( sup \N E \P + \{N E ); l A < C p ; e( sup [\Mf\P + \(M E ); l \r]) <C P ; 

\0<t<T / \0<t<T / 

E(\N E -N E \ P + \(N E ); 1 -(N E ); 1 \ P ) < C p \t-s\ p l 2 ; 

E(\Mf-M E \ p + \(M E ); 1 -(M E ); l \ p ) < C p \t-s\ p l 2 . 
Next, define 

1 



-E 



Z*ds\? Si _ 



Si- I 



Si-Si-1 

we are in the statement of Step 1 from which we deduce that 

£ / \Z E -ZZf\ 2 d S <C\n Q \. 
Therefore using again Lemma 3.4.2 of 11221 . page 71, we obtain 

HZ 

< IE 



/=i 



■Si— 1 



\Z?-Z?° \ 2 ds 

A Ji— 1 I 



/=i 
in 



Zf-Z*°> E \ 2 ds 



Si- I 
Si 



III n Si 

< £E / [\Z S -Z E \ 2 + \Z E -Zf^\ 2 }ds 



< E 



\Z,-Z E \ 2 ds 



+ C|7t |. 



(3.23) 
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Applying Proposition 2.4 we have 



limE 

e->0 



Jo 



\Z S -Z^\ 2 ds 



0. 



which, combined with ( 13.231 ), proves the theorem. 



□ 



4 Numerical scheme and rate of convergence 

In this section, we consider the BDSDE (12.21 ) in the special case ®(X) = h{Xj) where 
h G W 1 '°°(M. d ) such that h(0) is bounded by K. The goal of this section is to construct an 
approximation of the solution (X,Y,Z) by using the "step processes". Let recall ji : to < 

t\< < t„ = T the partition of [0, T] and \n\ = maxi< ( <„ | Af |, with Af = tj - f f _i. We 

set also A n Wi = W t , - W t ._, , A n B { = B t . - B tj { and for all < i < n define 

r,* = a(x t} , j<i)v?*, 

the discrete-time filtration. Let briefly review the Euler scheme for the forward diffusion X. 
Define Jt(f) = for t G Let X n be the solution of the following SDE: 

X* = x+ £b(n(s),X% s) )ds+ £a(n(s),X% s) )dW s , (4.1) 



and we define a "step process" X n as follows. 

x? = x% (t) , te[0,T\. 

The following estimate is well known (see e.g Kloeden and Platen, lfI31 ). 



(4.2) 



Proposition 4.1. Assume b and O satisfy the assumptions (HI) and (H3). Then there exists 
a constant C depending only on T and K, such that 



max E 

!</<« 



sup \X?-X t \ 2 + sup \X t -X t ._ 

0<t<T ?i-i<f<«i 



<c\%\. 



[21]. 



Moreover, we get the following estimate involving the step process X n due to Zhang in 



Proposition 4.2. Assume b and o satisfy the assumptions (HI) and (H3). Then there exists 
a constant C depending only on T and K, such that 

sup EflX/ 1 -X t \ 2 ] < C\%\; 

0<t<T 



E 



sup \X?-X t \ 

0<t<T 



< C|7t|log 



71 
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The backward component (Y,Z) will be approximated by the following numerical scheme: 
YZ = h(X?),Z«=0 



2^ = — EJL^A'W,] 



(4.3) 
(4.4) 



I* = Ef_j +f(t i - h Xl 1 ,Y* 1 ,Z«_ 1 )Al 
where E? [.] = E [. | and ?* = Y* + g(n , X% , F") A n B,. 

Remark 4.3. (?) Our approximation scheme differ from the one appearing in [6]. In- 
deed, actually we use the conditional expectation with respect the enlarge filtration 
J < i) V jF f s , which is necessary to extend Ito representation theorem for back- 
ward doubly SDE (see Pardoux and Peng, [20]). 

(ii) The backward component and the associated control (Y,Z), which solves the backward 

doubly SDE, can be expressed as a function of X and B, i.e. (Y t ,Z t ) = (u(t,B t ,X t ),v(t,B t ,X t 
for some deterministic functions u and v. Then, the conditional expectations, involved 
in the above discretization scheme, reduce to the regression of Y t K and Y t n (W tj — W tj _ , ) 
on the random variable {X^ ^B tj l ) . 

Next, for all < i < n, on can show that Y t n belongs to L 2 {Q.,f ti ), thus an obvious 
extension of Ito martingale representation theorem yields the existence of the (!F s ) se \ ti _ 1)ti \- 
jointly measurable and square integrable process Z n satisfying 



Y* = E[Y*\!F i K l ]+ f' ZfdW s . 

Therefore we define the following continuous version 

K-i " (' " ti-Mti-uX^ , r* , ,2* , ) - gfaX*, Y*){B t - B ti _ 



(4.5) 



y7t 



+ 



ZJdW, 



U-\ < t < U. 



(4.6) 



ti-i 



Note that the process Z K is given by the representation theorem, thus it is useful and even 
necessary to find a relationship with Z n define by (14.31) . We have 

Lemma 4.4. Assume b, o,f and g satisfy the assumptions (HI), (H2) and (H3) and let 

h 6 W l, °°(M. d ) such that h(0) is bounded by K. Then for all 1 <i <n, we have 

1 



Z % 



Zfds 



Proof. Let recall 



A K 7 K 



^^[[Y^ + gituX^Y^A^Bi) A n Wj] 



Then it follows from (031) that 



Zf = — E? , 
The result follows by Ito's isometry. 



A n Wi / ZfdW s 



□ 
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We also need the following which is the particular case of Theorem 3.2. 



Lemma 4.5. Assume b, o,f and g satisfy the assumptions (HI), (H2) and (H3) and let 

h G W L °°(R d ) such that h(0) is bounded by K. 
Let define, for each 1 < i < n, 



Z n 



Z s ds 



tt-i 



Then there exists a constant C depending only on T and K, such that 



E 



max sup 

l<i<nf i _ 1 <f< f[ 



i=l J tt-i 



<C\n\ 



(4.7) 



We are now ready to state our main result of this section, which provides the rate of 
convergence of the numerical scheme d4.3l )-( |4~4"l ). 



Theorem 4.6. Assume b, o,f and g satisfy the assumptions (HI), (H2) and (H3) and let 

h € W l, °°(M. d ) such that h(0) is bounded by K. Then there exists a constant C depending 
only on T and K, such that 



sup E\Y,-Y*\ 2 + E 

0<t<T 



\Z S -Z*\ 2 ds 



<c\%\. 



Proof. The proof follows the step of proof of Theorem 3.1 in Q so that we will only 
outline. In the sequel, C > will denote the generic constant independent of i and n and 
may vary line to line. For i £ {0, ...,n — 1}, we set 

5% = Y t - Y?, 5 % Z t = Z t -Zf, 8V(0 = f(t,X t ,Y t ,Z t )-f( ti ,X*,Y*,Z*) 
and &*g(t)=gM,Y t ,)-g(t i+1 ,Xl 1 ,Y* l ), t G [t h t i+l ). 

By Ito's formula, it follows from Lipschitz condition on /, g and h, together with the 
inequality ab < $a 2 + b 2 /fi that 



V t = E\§ n Y t \ 2 + E J'' +[ \5 K Z s \ 2 ds-\5% 



ti+i 



2E 



(8<Y s ,8*f(s))ds+ \5 K g(s)\ 2 ds 



ti+l 



< - £ ' +l E { \n\ 2 + \X S - X«\ 2 + \Y S - Y*\ 2 + \Z S - Zf} 2 } ds 
+ J"" 1 CE { |7i| 2 + \X S - Xf M | 2 + \Y S - Y£ +l I 2 } ds 
+$£' +l E\5 n Y s \ 2 ds, te[ti,t i+1 ). 



(4.8) 
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Proposition 3.1, Lemma 3.3 and Lemma 3.4 yield that 

E\X s -X£\ 2 + E\X s -X£ +l \ 2 <C\k\, 

E\Y S -Y*\ 2 < 2(E|F,-F f ,| 2 + E|8 7 %| 2 ) < C + E|8 7 %| 2 ) 
E\Y S - Y* +1 1 2 < 2 (E\Y S - \ 2 + E\8% +l | 2 ) < C (|7t| +E\8% +1 1 2 ) 

E\Z S - Zf} 2 < 2 (E\Z S - Z*| 2 + ^ Jl w E\5 K Z r \ 2 dr^ 
Plugging (|4~9|) into (|4~8|) . we get 

V t < ^f'^E^ + ^Ytf + ^-Zf^ds 

rti+i 



+cJ t ' +l E{\n\ + \S%J 2 }ds 



+^f t ' +l E\S K Z s \ 2 ds + $£' +1 E\8 n Y s \ 2 ds, 
from which and the definition of V t provide, for tj<t < , 

fti+i „ [-ti+\ 



E\b n Y t \ 2 + E\h n Z s \ 2 ds < pJ' +l E\8*Y s \ 2 ds+Ai 



where 



(l+C7i)E|5 ,t 7, +1 | 2 + 

+| ['' +[ E\d K Z s \ 2 ds. 
P Jti 



\n\ 2 + \n\E\Y t n \+ f' +l E\Z s -Zf. 
Jti 



Next, by little calculus used Gronwall's Lemma, we have 

E\8 n Y t \ 2 + J'' +l E\8 n Z s \ 2 ds < (1 + CP|jc|)A/; 
hence for t = ?,• and p sufficiently large than C, such that ^ < 1 , we obtain 

E|8 7t y f ,| 2 + (l-^)£ +1 E|5 Jt Z i | 2 ^ 

< (1 +C|7t|) |e|5% +1 I 2 + |7t| 2 + £ +i E[\Z S -2*\ 2 ]ds 

for small \n\. 

Iterating the last inequality, we get 

E\8%\ 2 + (\-^) f' +l E\h n Z s \ 2 ds 
P Ju 

< (l+C\n\) T ^\E\d*Y T \ 2 + \%\+j^ f' E[\Z s -Zf.J 2 ]ds 
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Moreover, it follows from Lemma 3.4, Lipschitz condition on g and Proposition 3. 1 that 

C fti+i. 



E\d%\ 2 + (1 - [ ' +> E\S K Z s \ 2 ds 
p Jti 

< (l+C\%\) T/% {E\8 K Y T \ 2 + \%\+C\%\} <C\n\ 

for small |tc|. 

On the other hand summing up inequality (14.1 II) with t = t{, we get 



i-|(i+cpN) 
c, 



[ T E\h n Z K \ 2 ds 
Jo 



< (1 +Cp|jc|) J|jc| + (1 +Cp|7t|)(l +C|7t|)lE|6 7t Fr| 2 



+ 
+ 



(i+cpN)^N-i 



ElS^ol 2 



C, 



i n-l 



(i+cpw)((i+cN)+p|7u|)-i 



£E|8 n y fi |2 



i=l 



+(i+cpN)^£ 

P i=0 Jt i 



E\Z s -Z*\ 2 ds. 

Therefore, by inequality (14.121 ) and Lemma 3.4 one derives that 

[ T E\h n Z s \ 2 ds<C\%\ 
Jo 



(4.12) 



and then 

sup |8*y,| 2 <C\n\. 

0<t<T 

□ 

To end this section, let give the following bound on K^'s which will be used in the 
approximating of discrete conditional expectation Ef, for all < i < n — 1. 

Lemma 4.7. Assume b, o,f and g satisfy the assumptions (HI), (H2) and (H3) and let 

h 6 W l, °°(W 1 ) such that h(0) is bounded by K. For all < i < n — 1, define the sequences 
of random variables by backward induction 

a* = 2C, p„ = C, 

a? = (l-C|7i|)- 1 (l+C 2 |7i|) 1/2 {(l + 2C|7i|)[(l+C|A It B,- +1 |)Pf +1 +C|A It B ; - +1 |]+C|7i|} 
Pf = (1 -C\n\)- l (\ +C 2 \%\) l l 2 {(1 +C| A*B i+l \)a? +1 +6C 2 \n\(l + 2C\ A«B i+1 \) + 3C\n\} 

Then, for all < i < n; 



\Y t f\ < a? + $?\X* 12 



(4.13) 
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< (1 + 2C|7U|) {(l+C\ A n B i+ i\)$f +l +C\ A*B i+1 \} \X*\ 2 

+(1 +C\ A n B i+l \)a? +l +6C 2 |7i|(l + 2C| A % B i+l \)labelbound2 (4.14) 

rati I (y* +*M«, nf) A 71 fl f ) A 71 w i+l I 

< v^f(l+2CN){(l +C| A"B i+l \)&? +1 +C\ A«B i+1 \} \X«\ 2 
+y/\n\(l+C\ A 7C B ;+1 |)af +1 +6C 2 |7i|(l + 2C| A K B i+l \). (4.15) 

Moreover, 

limsup max (ocf + Pf ) < °°, a.s. 

0<i<n 

Proof. First, since h is C-Lipschitz, h(0) is bounded by C, 

\Y tn \ = \h(X?)\ <C(\X?\ + l)<2C + C\X?\ 2 = a n + p„|X*| 2 . 
Next, we assume that 

|F, +1 |<af +1 +Pf +1 |^J 2 , (4.16) 

for some fixed < i < n — 1. Then by the definition of Y n in (14.41) . there exists a im- 
measurable random variable Q such that 

(1-C|jt|)|y»| < E?[(y£ 1+ ^ +1 ^^^ 

< (Ef|7 J f +1 +g(f ( - +1) X^,y f f + JA w B i+ i| 2 )V2(Ef| + C-A^- +1 | 2 ) 1 /2 

+C|7i|(3 + |^| 2 ). (4.17) 

It show in [ 6 ] that 

Efll + ^A^+il 2 < 1+C 2 |7i|. 

This provide from (14.171) 

(l-C|7t|)|F*| < (l + C 2 |; t |) 1 /2(Ef|F^+g( fi+ i,X^ ) ^ 1 )A' t B i+ i| 2 ) 1 /2 

+C|7i|(3 + |X f 7C | 2 ). (4.18) 
But it follows from the Lipschitz property of g that, 

mK.+g^uX^JA-B^l 2 ) 1 / 2 

< (1 + C\ A n B i+{ |)(E? |y f f I 2 ) 1 / 2 +C| A** m |(Ef |X* , I 2 ) 1 / 2 



< 



(l+C|A It B ; - +1 |){af +1 +Pf +1 [(l + 2C|7 I |)|^| 2 + 6^ 2 |7 I |]} 
+C|A 7t B,- +1 |[(l+2C|7t|)|X^| 2 + 6^ 2 |7t|] 
(1 + 2C|ju|){(1 +C\ A^,- +1 |)Pf +1 +C| A"B,+i|} \X*\ 2 
+(1+C| A 7C B m |)a? +1 +6C 2 |7i|(l+2C| A n fl, + i|). 
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Finally (14.181) becomes 

\Y t *\ < (i-cItiD-^i+c 2 !^) 1 / 2 

x {(1 +2CW)[(1 +C| A*B i+1 \)ft +l +C\ A K B i+ i\] +C\n\} \X*\ 2 

+ (1-C|7I|)- 1 (1+C 2 |7I|) 1 / 2 

x {(1 +C\ A n B i+1 \)a* +1 + 6C 2 |ti|(1 + 2C\ A K B i+ i\) + 3C|tc| } 
= a? + (3?|X*| 2 . 

□ 

5 Rate of convergence of the regression approximation 

In this section, we try to give some ideas for method of simulating numerical scheme de- 
rived in the above section. It well know that the process X % defined by d4.il ) is simulated by 
the classical Monte-Carlo method. We are reduced to simulate the process (Y n ,Z K ) defined 
in (14.31 ) and ( I4.4I ). In practice, the main tool to define of an approximation of Y n , and then of 
Z 71 , is to replace the conditional expectation Ef by its estimator Ef in the backward scheme 
(14.31 ) and (14.41 ). We first establish the following bound on the F^'s which help us to derive 
this simulation. 

For the regression approximation, we consider {&f}Q<i<n, {5Lf}o<i<« an d {j?}o<i<n 
defined by: 

a? + Pf|X,f| 2 

(1 + 2C|jt|) {(1 +C\ A*B i+l \)ff +1 +C\ A % B i+ i\} \X«\ 2 
+ {l+C\A 7Z B i+l \)af +l +6C 2 \n\(\ + 2C\A K B i+l \) 
v/R(l + 2C|7t|) {(1 +C\ A*B i+1 \)tf +1 +C\ A*B i+1 \} \X*\ 2 
+^(l+C|A 7t B i - +1 |)af +1 +6C 2 |7i|(l+2C|A 7t B ; - +1 |). 

Therefore thanks to Lemma 4.7, we have 

- af(X* A*B i+1 ) < Y* < P?(X*, A K B i+l ) (5.1) 

A*B i+1 ) < Ef [F*J < n?(X*,A*B i+1 ) (5.2) 

-J?(X*, A*B i+i ) < E[F* ( A 71 W i+ i] < tf(X*, A K B i+i ) . (5.3) 
Next, for a M-valued random variable we define 

if® = -2?(X*,A*B i+1 )\/1iA*?{X«,A*B i+1 ) 

T*"(£) = -^(X^A n B i+1 )W^A3if(X^A K B i+l ) 

Tf® = -Jp(X^,A K B i+ i)y^/\jP{X^,A K B M ). 



v? = 
= 

jf = 
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Given an approximation Ef of Ef, we are ready to get the process (Y n ,Z n ) defined by 
following backward induction scheme: 



I t„ 


= g&¥), 


(5.4) 


7* 


= ^Ef_ j [ (t£ + g(tj , X* , tg) A n Bi) A"W f ] 


(5.5) 


Vi 


= E^ + gM* A'BJ + /(* ! -_ 1 ,X£ l ,F f ?_ i ,Z£_ 1 ) Af 


(5.6) 


v n 


= (Y n ) 


(5.7) 



for all 1 < i < n. 

Remark 5.1. Using the above notation and replace X K by U n = (X K ,B n ), on can state anal- 
ogous to examples 4.1 and 4.2 appearing in |6j. 

To end this section, we derive the following LP estimate of the error Y % — Y n in terms 
of the regression errors Ef — Ef. 

Theorem 5.2. Let p>\ be given, and (P n be a sequence defined above. Then, there is a 
constant C depending only on T,K and p such that 

+ ||^y-Ey)[(^ 1 +^ 1> l* 1> ^JA"ll /fl )A*W H .l]|| W } 

Proof. For < i < n — 1 be fixed, with a similarly calculus as one use in proof of Theorem 
4.1 (6l, we have 

(1 - Cn)\Y t f -f£\< N + (1 + C| A n B i+l \)(E\Y t * - y^| p ) 1/p (Ef|l + Q A" W i+l \ 2k fl 2 \ 

(5.8) 

where k is an arbitrary integer greater than the conjugate of p and 

+^ +1 {f(t h xiY t *,(A? +l )- l m?^ 

-f(t h XlYZ,(A? +1 )- l 2j[(YZ +l +g(t i+h XlJZ l )A*B i+1 )A*W M })}. 

Since 

\\£i\\ LI > < Tli = C( || (E,- - E,-) + g(t i+l ,X* , , % , ) A^B i+l ]\\ u 
+||(E ( --EO[(^ 1 +g(f/ + i,^ l! f«JA w B i - +1 )A' c ^ +1 ]||„, 

and 

(Ef|i + ^A 7t w;- +1 | 2 *) 1 / 2 *<(i+cN) ) 
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we get from (15.81) 



(l-CK)\\Y t f-Y t f\\^<m + (l + C\K\) 





LP 



(5.9) 



and thus the result follows as in [6]. 



□ 



Remark 5.3. With the foregoing results, it is possible with some not very difficult adjust- 
ment to obtain similar results as those obtained by Bouchard and Touzi (see Section 5 and 
6 of El). 
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